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I INTRODUCTION 


Recent  results  have  demonstrated  that  the  use  of  recei ve-array 
beamforming  systems  which  are  designed  to  be  “optimal"  in  the  sense  of 
rejecting  noise  and  interference  offer  superior  performance  over  con- 
ventional beamforming  methods  in  HF  applications.  The  purpose  of  this 
report  is  to  offer  a quantitative  analysis  of  the  effects  of  random 
fluctuations  in  the  received  waveforms  on  the  overall  performance  of 
optimal  beamformers.  These  fluctuations  are  the  result  of  signal 
distortion  imposed  by  ionospheric  reflection  as  well  as  the  presence 
of  system  noise  in  the  recei ve-array  electronics.  In  this  report,  we 
assume  that  signal  distortion  affects  both  gain  and  phase  of  each 
spatially-propagated  waveform  received  by  the  elements.  Thus,  the 
effects  of  gain  and  phase  errors  on  both  the  desired  or  target  signal 
as  well  as  directional  interferers  are  considered.  In  the  case  of 
phase  errors  on  the  desired  signal,  we  consider  two  separate  cases: 

1)  the  zero-mean  case  in  which  the  average  direction  of  arrival  of  the 
target  signal  coincides  with  the  pointing  direction  of  the  array  and 

2)  non-zero  mean  cases  in  which  the  mean  arrival  direction  is  other 
than  expected.  In  the  latter  case,  the  desired  signal  is  treated  as 
an  interfering  signal  by  the  processor  and  may  suffer  some  amplitude 
degradation  at  the  beamformed  output.  Quantitative  predictions  of  this 
effect  as  a function  of  the  offset  angle  are  presented  in  section  IV  of 
this  report.  Comparisons  with  conventionally-weighted  beamforming 
systems  are  also  presented  for  all  cases  studied. 

Previous  studies  relating  to  the  effects  of  random  gain  and  phase 
errors  have  largely  been  restricted  to  the  case  of  conventional  beam- 
forming systems.  (See  Ruze  [1],  Rondinelli  [2],  and  Green,  et.  al.  [3]) 
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Carver,  et.al.,  [4]  have  studied  the  effects  of  beam-pointing  errors 
in  planar  phased  arrays  and  recent  work  by  McDonough  [5]  has  shown 
the  effects  of  random  input  errors  on  the  directivity  of  both  con- 
ventional and  linearly-constrained  optimal  beamfortners.  The  work 
presented  in  this  report  augments  that  of  McDonough  in  that  two 
different  optimal  structures  are  analyzed  herein  and  specific  cal- 
culations of  noise  rejection  and  beam  patterns  are  presented. 

In  section  1 1 , we  describe  the  conventional  and  optimal  beamform- 
ing structures  of  interest  in  this  work.  The  latter  are  termed  noise- 
cancelling beamformers  and  utilize  auxiliary  beamforming  subsystems 
designed  to  cancel  undesired  interference  which  is  present  in  the  con- 
ventional array  output.  The  specific  model  used  to  compute  the  effects 
of  random  errors  is  described  in  detail  in  Section  111.  Briefly,  we 
model  each  spatially  propagated  signal  incident  on  the  array  as  having 
been  transmitted  to  the  array  by  a random  channel.  The  characteristics 
of  this  channel  then  determine  the  nature  of  the  phase  and  gain  errors 
observed  at  the  element  outputs.  Simulation  results  based  on  this 
model  are  presented  in  Section  IV  and  Section  V contains  summary  and 
conclusions  drawn  from  the  study. 


II  BEAMFORMING  STRUCTURES 


The  general  form  of  the  receiving  array  stem  of  interest  in  this 
work  is  shown  in  Fig.  1.  Each  array  element  contains  a separate 
receiver  which  converts  the  HF  element  signals  to  baseband.  In  practice, 
subarrays  may  serve  as  elements  for  each  individual  receiver  [6,7]  and 
the  individual  receivers  are  gain  and  phase  matched  across  the  band- 
width of  the  transmitted  signal.  The  baseband  outputs  are  then  simul- 
taneously sampled  and  digitized  prior  to  beamforming.  A sampling  rate 
greater  than  twice  the  receiver  bandwidth  is  used  to  prevent  aliasing 
and  all  beamforming  is  carried  out  using  a real-time  digital  processor. 

We  assume  that  the  receivers  shown  in  Fig.  1 contain  steering  delays 
which  pre-point  the  array  system  to  the  assumed  direction  of  the  desired 
signal.  Thus,  in  an  ideal  system,  if  x.(k}  is  used  to  represent  the 
k-th  sample  from  the  i-th  receiver  output, 

X.  (k)  -=  s(k)  + n.  (k)  (1) 

where  s(k)  is  the  desired  signal  and  n.(k)  is  the  total  of  all 
noise  and  interference  observed  at  the  i*"^  element.  Of  course,  in 
practical  systems,  the  presence  of  system  and  propagation  errors  will 
cause  the  s(k)  term  in  Eq.  (1)  to  be  a function  of  element  number  i. 

Conventional  beamforming  consists  of  gain  weighting  the  receiver 
samples  x.(k)  to  produce  a beamformed  output  given  by 

y.(k)  = I 

C • t ^ I I 

1 = 1 

where  W . is  the  i-th  conventional  weight. 

Optimal  beamforming  structures  employ  a beamformer  which  utilizes 
delayed  samples  from  each  receiver.  These  processors  have  been  proposed 
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Element 


Figure  1.  Digital  Antenna  Array  System  Model 


Figure  2.  Software-constrained  noise-cancelling  optimal  array  structure 
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in  a variety  of  forms  and  have  been  widely  discussed  in  the  literature 
[6-11].  The  specific  structures  of  interest  in  this  work  have  been 
termed  noise-cancelling  or  side-lobe-cancelling  array  structures  [12,13]. 
Within  this  class  we  consider  two  distinct  implementations  termed 
software-constrained  (Fig.  2)  and  hardware-constrained  (Fig.  3). 

Software-Constrained  Optimal  Beamformer 

In  the  software-constrained  system,  a conventional  beam  output 

y (k)  is  formed  in  a manner  identical  to  that  shown  in  Eq.  (3)  except 
c 

for  the  presence  of  a bulk  time-delay  of  1/2  samples.  Thus, 


y (k)  = y w .x.(k-L/2) 

c'  Cl  i' 


(3) 


The  reason  for  the  delay  will  be  explained  shortly.  A noise-cancelling 
auxiliary  beam  y (k)  is  formed  from  a weighted  sum  of  delayed  signals 

d 

from  each  of  the  receivers.  A total  of  L samples  from  each  element 

are  combined  using  the  multichannel  digital  filter  shown  in  Fig.  4. 

in  matrix  notation,  we  denote  the  vector  of  KL  auxiliary  weights 

as 

-a 


S 


S 

;I,L 

w^ 

^k,l 

'S 

^k,L 


W'  = 


(4) 


These  weights  multiply  a vector  of  delayed  receiver  signals  X(k)  defined 
by 
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Figure  3. 


Hardware  constrained  noise-cancelling 
optimal  array  structure 
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y,(k)  = X^(k)W^  (6) 

a 


In  noise-cancelling  array  structures,  an  optimal  weight  set 
is  found  such  that  the  output  noise  and  interference  power  is  minimized 
subject  to  the  constraint  that  the  desired  signal  passes  undistorted 
through  the  processor.  Since  the  desired  signal  is  presumed  to  be 
in-phase  at  each  receiver  output — Eq.(l) — this  can  be  accomplished  by 
constraining  each  column  of  weights  in  the  processor  shown  in  Fig.  k 
to  sum  to  zero.  Thus  we  require 


Multichannel  Digital  Filter 


in  this  expression,  ^ 9 have  been  used  to  denote  column  vectors 

of  ones  and  zeros,  respectively. 

Once  the  weighting  coefficients  have  been  constrained  by  Eq.  (8) — 
termed  here  the  software  constraint'- the  ideal  signal  cannot  appear  in 
the  auxiliary  beamformed  output.  Minimum  output  noise  and  interference 
is  then  achieved  by  finding  such  that  the  difference  beam  y^(k) 

in  Fig.  2 has  minimum  power  subject  to  the  constraint  in  Eq.  (8). 

Frost  [1 4]  has  given  an  excellent  discussion  of  this  problem  and  the 
resulting  solution.  As  discussed  in  [6,7]>  inserting  a delay  of  L/2  in 
the  conventional  beam  places  the  phase  reference  of  the  conventional 
beamformed  signal  in  the  center  of  the  auxiliary  processor.  (Note  that 
if  L is  odd,  L/2  is  replaced  by  (L-l)/2).  This  phase  shift  has  been 
shown  experimentally  to  offer  the  best  overall  performance  in  practical 
implementations  of  noise-cancelling  arrays. 

In  order  to  obtain  a convenient  expression  for  the  form  of  the 
optimal  weights  W^,  we  defined  a KL-dimensional  conventional  beam- 
forming vector  as 

g 

0 


Wc  = 


w 


cl 


w 


ck 


(10) 


The  optimal  weights  are  then  given  by  [11,  14] 


- w - r*'c(c^r“'c  )'*C^W 
-a  -c  "x  = * »x  =x  = -c 


(ID 
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where 


R is  the  KL  x KL  matrix  of  space/time  correlations  observed 
=x 

at  the  receiver  outputs.  Specifically, 

R^  ^ E [X(k)X^(k)l  (12) 

where  E[*]  denotes  expectation. 

In  practical  applications  of  this  beamformer  such  as  those  dis- 
cussed by  Kelly,  e^.aj_.  [11],  the  correlation  matrix  is  computed  using 
of  observed  receiver  output.  For  example,  the  correl- 
the  signal  from  element  i delayed  by  H units  and  that 
j delayed  by  m units  is  approximated  by 
, T 

r.  .()l,m).  ^ ^ X.  (k-£)x.  (k-m)  (13) 

I iJ  ' J 

° k-1 

^ is  the  observation  interval.  For  the  purposes  of  the 

present  study,  we  assume  that  averages  of  this  or  an  equivalent  type 

have  been  carried  out  with  sufficient  data  to  accurately  estimate  the 

S 

elements  of  R . This  matrix  is  in  turn  used  to  compute  W using 
Eq.  (11) 

Hardware-Constrained  Optimal  Beamformer 

Figure  3 shows  an  alternate  implementation  of  the  s idelobe-cancel 1 ing 
optimal  array  structure.  The  only  difference  between  this  and  the  software- 
constrained  beamformer  in  Fig.  2 is  the  use  of  a spatial  preprocessing 
matrix  which  converts  the  vector  of  received  signals  X(k)  into  a 
modified  vector  X'(k).  The  auxiliary  beam  y (k)  is  formed  as  a 
weighted  sum  of  the  elements  of  X'(k)  using  a multichannel  filter 


similar 

to  that  shown 

in  Fig.  4. 

The  purpose  is  to  block  the 

des i red 

signal  s(k) 

from  X' (k) 

where 

time  averages 
at  ion  between 
from  element 


where  T 
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X'(k)-W^X(k)  (lA) 

is  a K'L  dimensional  vector  of  modified  receiver  signals  (K*  < K) . 

Inspection  of  Eq's  (1),  (5),  and  (14)  shows  that  s(k)  will  be 

successfully  blocked  if  is  constructed  as  a diagonal  matrix  of 

identical  K'  x K-dimensional  submatrices  w , 

=s 


w 0 0 ...  0 

ms  ■ • m 

0 w 0 

m ms  m 

W 

■S 

where  0 is  used  to  represent  a K x K matrix  of  zeros.  In  effect, 
implementing  the  preprocessor  in  the  form  shown  in  Eq.  (15)  is  equivalent 
to  applying  a spatial  prefilter  w^  to  the  set  of  K received  array 
signals  {x.(k)},  and  s(k)  is  blocked  provided  that  the  rows  of  w^ 
sum  to  zero. 

Two  examples  of  spatial  preprocessors  for  the  case  where  K ■ 8 
are  shown  in  Fig.  5.  In  Fig.  5a,  K'  ■ 4 and  w^  is  given  by 


Alternately,  in  Fig.  5b,  K'  - 7 and 
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(17) 


Other  forms  of  spatial  preprocessors  which  successfully  block  the 
desired  signal  are  possible  and  discussed  further  in  the  sections 
following. 

Given  that  X'(k)  does  not  contain  s(k) — i.e.  the  desired 

signal  has  been  blocked  by  the  hardware  preprocessoi — minimum  output 

noise  and  interference  is  achieved  by  finding  W such  that  the 

difference  beam  yQ(k)  in  Fig.  3 has  minimum  power.  The  equation  for 

the  auxiliary  beam  y (k)  in  matrix  notation  is 

d 

y“(k)  - x^(k)wV  (18) 

o * *5  a 


and  the  resulting  optimal  weights 


are  given  by  { 


-a 


W R W 


*S“X-C 


(19) 


A complete  derivation  of  this  result  is  presented  in  Appendix  A.  It 
should  be  noted  that  the  correlation  matrix  R^  which  appears  in  (19) 
is  identical  to  that  defined  by  Eq's  (12)  and  (13)  for  the  software- 
constrained  beamformer. 

Discussion 

It  has  been  demonstrated  [II,  15]  that  the  optimal  software  and 
hardware  auxiliary  beamformers  shown  in  Figs.  2 and  3 and  defined  by 
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Eq's  (11)  and  (19),  respectively,  are  identical  provided  that  the 
spatial  matrix  is  of  rank  K-1.  Equivalently,  K‘  * K-1  and  the 

rows  of  must  be  linearly  independent.  One  may  interpret  this 

hardware  constraint  preprocessor  as  a set  of  K'  beamformers,  each 
having  a null  in  the  direction  of  the  desired  signals.  Methods  for 
encoding  additional  constraints  into  the  preprocessor  are  discussed 
in  Reference  [12]  and  [13l- 
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PERTURBATIONS 

The  formulation  presented  in  the  previous  section  for  optimal 
array  structures  assumes  that  the  desired  signal  is  a plane  wave 
incident  on  the  receive  array  from  a known  direction.  Under  these  con- 
ditions, optimal  array  processing  can  be  shown  to  provide  significant 
improvement  in  the  rejection  of  undesired  directional  interference 
when  the  interference  consists  of  plane  waves  propagating  uniformly 
across  the  array.  However,  when  the  target  signal  comes  from  other 
than  steered  direction,  the  optimal  processor  may  regard  the  signal 
as  interference  and  present  some  degree  of  rejection  of  the  desired 
signal.  In  addition,  some  degradation  of  performance  may  result  due 
to  amplitude  and  phase  fluctuations  of  all  spatially  propagated 
signals.  The  purpose  of  this  section  is  to  present  a model  which  can 
be  used  to  quantitatively  determine  the  effects  of  these  errors  on 
the  overall  performance  of  both  forms  of  the  optimal  processor. 

In  the  following,  we  will  assume  a stationary  environment  and 
model  the  amplitude  and  phase  variations  as  zero-mean  random  processes. 
The  effect  of  an  error  in  the  assumed  direction  of  arrival  of  the 
desired  signal  is  included  by  assuming  a non-zero  mean  for  the  phase 
errors.  For  simplicity,  this  mean  corresponds  to  plane-wave  propaga- 
tion from  other  than  the  desired  direction  and  performance  degradation 
is  measured  as  a function  of  angular  offset,  or  tilt,  of  the  desired 
signal.  In  addition,  it  is  further  assumed  that  the  antenna  array 
consists  of  K elements  distributed  in  a line  at  equispaced  intervals. 
The  direction  of  arrival  of  the  desired  signal  is  perpendicular  to 
the  array  and  is  termed  the  bores ight  direction. 
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Model  and  Assumptions 


To  simplify  the  following  presentation,  the  case  of  a single 
spatially-propagating  desired  signal  and  a single  directional  inter- 
ference is  first  considered.  The  model  formulation  is  then  extended 
to  the  case  of  multiple  signals  and  interferences. 

Figure  6 shows  the  signal,  noise  and  interference  model  for  the 
simple  case.  Two  spatially  coherent  signals  are  present:  a desired 
target  signal,  s(t),  and  a directional  interference,  n(t).  They 
are  modeled  as  zero-mean,  random  and  uncorrelated  processes,  and  are 
assumed  to  be  sinusoidal  with  random  phases.  Two  independent,  identi- 
cally distributed  channels  denoted  by  gQ  and  gj  in  Figure  6 are 
used  to  model  the  effects  of  random  amplitude  and  phase  fluctuations. 
Thus,  the  output  signals  s(t),  n(t)  contain  random  errors  and  pro- 
pagate toward  the  receiving  array  in  the  form  of  distorted  plane  waves. 
To  illustrate  the  channel  model  consider  the  target  signal  shown  in 
Figure  6.  For  a complex  signal  representation,  we  have 


s(t) 


a e 
s 


. (u)  t+<l)  ) 

J O ^s 


(20) 


where  and  are  the  amplitude  and  demodulated  frequency,  res- 

pectively, of  the  desired  signal.  The  phase  term  (j)^  is  modelled  as 
a uniform  random  variable  over  a 2n  interval.  The  power  content 
of  this  signal  is  given  by 

E[|s(t)|^]  = of  (21) 


The  random  channel  adds  both  gain  and  phase  modulation  to  this 
signal  and  produces  an  output  s(t)  given  by 
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(22) 


pr  ' f 


s(t)  = c^[l+a(t)]a^e 


Both  amplitude  and  phase  variation  terms  a(t)  and 
as  independent,  zero-mean,  gaussian  processes  with 
respectively.  The  scalar  term  c^  is  included  to 
of  the  input-out  power  of  the  channel  and  is  given 


0 (t)  are  model  1 ed 

2 2 

variances  a and  a. 

a 0 

preserve  invariance 
by 


c 

s 


(23) 


The  time  behavior  of  a(t)  is  modelled  as  white — i.e.  wide  band — 
and  0(t)  is  assumed  to  be  constant  over  the  time  interval  during 
which  correlations  are  measured.  Each  array  element  is  presumed  to 
receive  an  s(t)  which  results  from  a separate  channel  which  is 
independent  of  the  other  array  elements.  In  effect,  this  model  assumes 
that  amplitude  variation  are  uncorrelated  spatially  and  rapidly  time 
varying  but  that  the  phase  variations  are  spatially  decorrelated  and 
nearly  constant  in  time  at  each  element.  Assumptions  identical  to 
these  have  been  used  by  Lee  and  Waterman  [16]  to  accurately  model 
millimeter  wave  propagation  through  the  atmosphere. 

Thus,  after  propagation  through  the  i-th  channel,  the  sampled 
signal  Sj(k)  observed  at  the  output  of  the  i-th  receiver  is  given  by 


s.  (k) 


c^[l+a. (k)] 


j[aj  k+(l)  -0  (k)] 

o e 
s 


The  generalized  spatial/tempral  correlation  element  defined  by  Eq. 
is  then  modelled  for  the  signal  only  and  (assuming  i!i<j)  as 


(24) 

(13) 


i 


1 6a 


r!  = E[s.  (k-&)Sj  (k-m)]  = E [s . (k-£)  ] E [sj' (k-m) ] 


j(u)  k-wJl+(l>  ) -j0.(k-il) 

= c^o^e  ° ° " E[e  ' ] 


-j(w  k-co-m+(J)^)  j0  (k-m) 

. c a e ° ° ^ E[e  J ] 

S 5 


(25) 


-j0j(k-Jl)  j9.(k-m) 

But  the  two  terms  E [e  ]E[e  ] are  recognized  as  the 


characteristic  functions  of  gaussian  random  variables  [17]  and  each 

.2 


-0 


0/2 


each  have  value  e • The  term  in  (25)  then  reduces  to 


s , 2 2 -<^0 

r.  . = c n e e , 

\ yj  5>  b 


(26) 


When  i=j  is  readily  shown  that 


s , 2 2 

ri  (£.01)  = c^a^e 


, £ ?^  m 


= o 


£ = m 


(27) 


For  the  case  of  real  sinusoids,  the  correlation  expressions  in  (26)  and 
(27)  are  replaced  by  their  equivalent  real  parts.  Further,  by  defin- 


ing a gain  reduction  factor  as 


“s  " , . 2 


1+0 


(28) 


and  using  (23),  the  resulting  spatial/temporal  correlation  functions 


for  a real  signal  may  be  expressed  as 

2 


r^  .(£,m)  = 0^(l-ot  )e  °cos[co^(m-£)] 
I » j s a u 


rf  ;(£,m)  - a^(l-a  )cos  [<i)„(m-£)] 
* » * ^ a u 


i / j 

£ m 


(29) 


£ = m 


(30) 


1 


18 


Discussion 

The  development  presented  above  in  Eq's  (24)-(30)  has  shown  how 
correlation  values  for  the  desired  signal  may  be  computed  from  the 
random  channel  model  shown  in  Fig.  6 and  described  by  Eq's  (22)  and 
(23).  For  a directional  interference,  analagous  calculations  are 
carried  out.  The  only  difference  in  this  case  is  that  the  random 
phase  term  6.(k)  in  (2k)  is  modeled  as  having  a non-zero  mean  6. 
corresponding  to  the  average  direction  of  arrival  of  the  interfering 

wavefront.  The  term  E [e  ] in  (25)  then  becomes  [17]  i 

j6.  -^0/2 

e e , and  Q.  represents  the  average  phase  delay  to  the  i-th 
array  element.  Spatial/temporal  correlation  values  for  a directional 
interferer  similar  to  those  in  Eq's  (29)  and  (30)  are  then  given  by 


2 

'0l„ 


j(l,m) 

2 "“9  - - 

= 0 (1-a  )e  cos  [w„(m-Jl) +0.-0.] 

n a u j 1 

. i ^ j 

(31) 

j (Jl,m) 

2 

= Op cos  [ajQ(m-Jt)] 

, 1 ^ n 

2 

, i!,  = m. 

(32) 

The  random  channel  model  is  readily  extended  to  the  case  of 
multiple  directional  interferers  and/or  multiple  target  signals  by 
modeling  these  sources  as  statistically  independent.  In  this  case, 
the  net  correlation  function  is  the  sum  of  appropriate  signal  and 
noise  terms  similar  to  those  given  above. 

Inspection  of  the  correlation  terms  in  Eq's  (25)  and  (31)  shows 
that  the  effect  of  the  gaussian  amplitude  and  phase  noise  in  the 
system  is  to  reduce  the  amount  of  inter  element  correlation.  For 
the  amplitude  model,  this  reduction  term  is  (1-a  ) and  for  phase 
< 

noise  it  becomes  e . The  latter  is  a direct  consequence  of  the 
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gaussian  assumption  and  changes  as  other  statistical  models  are  used. 

For  example,  if  a uniform  probability  density  on  phase  is  used  with 

2 

zero  mean  and  variance  a , the  term  is  given  by  a sine  function 
sin(i/30Q)//30Q. 

Because  the  model  reduces  the  measured  spatial  correlation  values, 
one  would  expect  a reduction  in  the  amount  of  improvement  offered  by 
optimal  arrays  to  be  more  comparable  with  that  provided  by  conventional 
beamformers.  In  the  limit  of  ::omp1ete  spatial  decorrelation  it  is 
readily  shown  that  optimal  and  conventional  processors  have  identical 
performance.  Examples  presented  in  the  section  following  illustrate 
this  degradation  for  the  case  of  a simple  line  array  of  eight  equally- 
spaced  elements. 


20 


IV  SIMULATION  RESULTS 


In  this  section,  we  present  some  simulation  which  illustrate  the 
sensitivity  of  a line  array  to  random  gain  and  phase  errors.  The 
simulations  were  carried  out  using  the  CDC  6400  facility  at  the 
University  of  Colorado  Boulder  campus.  A receiving  array  of  eight 
elements  was  chosen  and  a raised-cosine  window  was  selected  for  the 
conventional  beamformer.  The  optimal  beamformers  are  characterized 
in  terms  of  number  of  spatial  and  temporal  degrees  of  freedom  (OOF)  in 
the  auxiliary  beamforming  path.  Thus,  a four  element  auxiliary  channel 
with  four  taps  per  delay  line  is  indicated  by  4 x 4 DOF. 

Input  Signal  and  Noise  Fields 

Two  cases  were  considered.  Case  (i)  consisted  of  a single  target 
signal  and  one  directional  interference,  and  case  (ii)  consisted  of 
multiple  target  signals  and  interferences.  In  both  cases,  uncorrelated 
receiver  noises  were  added  to  the  inputs  and  the  total  power  input 
normalized  to  unity  at  each  element.  The  angle  of  arrival  for  all 
directional  sources  was  measured  in  terms  of  a normalized  spatial  angle 
0 chosen  in  the  range  of  -tt  to  it  . An  angular  increment  of  Zn 
radians  corresponds  to  the  spatial  fold-over  angle  of  the  linear 
array.  This  angle  is  generally  referred  to  as  the  gratinglobe  separ- 
ation and  the  actual  angle  in  a practical  array  depends  on  the  element 
separation  and  HF  carrier  frequency  [6].  The  received  frequency  f 
used  for  the  simulations  is  presented  on  a normalized  scale  of  0 to 
0.6  with  0.6  representing  one-half  of  the  sampling  frequency. 

Table  I summarizes  the  details  of  the  input  waveforms.  The 
signal -to- interference  ratio  (SIR)  and  the  signal-to-total-noise 
ratio  (STR)  observed  at  each  element  were  defined  as  follows: 
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TABLE  1 


V 


SUMMARY  OF  INPUT  WAVEFORM  CHARACTERISTICS 
(a)  Single  Target  Signal  Example 


Target  Signal 

1 nterference 

Receiver  Noise 

Power 

f 

6 

L_ 

Power 

■ 

6 

Power 

0.02 

.25 

D 

CD 

00 

.20 

0.4 

1 

0.1 

(SIR) - -l6.43dB 


(b)  Multiple  Target  Signal  Example 


Target  Signal 

1 nterference 

Receiver  Noise 

Power 

f 

0 

Power 

f 

0 

Power 

0.01 

0.225 

0 

mm 

0.2 

0.1 

0.01 

0.25 

0 

H 

0.3 

0.01 

0.275 

0 

B 

0.4 

(SIR) . = -14.77  dB 
1 n 
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(S,R),„ 


total  input  power  of  all  target  signals 

total  input  power  of  all  directional  interference 


(STR)  . 
in 


total  input  power  of  all  target  signals 

total  input  power  of  directional  interference  and  receiver  noise 


An  evaluation  of  the  performance  observed  with  these  examples  was 

performed  and  is  summarized  below.  As  discussed  in  Section  II  above, 

the  software  and  hardware  constrained  beamformers  are  identical  when 

the  spatial  preprocessing  matrix  w (see  Fig.  3)  produces  K-1 

=s 

independent  outputs.  In  the  comparisons  which  follow,  then,  any 
optimal  beamforming  structure  which  is  characterized  as  having  7 x L 
OOF  applies  in  both  optimal  structures.  If,  on  the  other  hand,  K'  x L 
OOF  are  used  and  K'  < 7 then  the  results  apply  only  to  the  hardware 
constrained  case.  To  quantitatively  measure  the  performance  of  the 
array  structures  of  interest  in  the  presence  of  random  gain  and  phase 
errors,  four  criteria  were  selected.  These  were: 


A.  Main-Lobe  Gain  Reduction:  The  gain  reduction  in  the  main  lobe 
direction  was  computed  as  the  ratio  of  output  target  signal  power 
under  random  perturbations  to  that  under  no  perturbation.  Figure  7 
illustrates  this  parameter  as  a function  of  both  the  amplitude  noise 
variance  reduction  factor  oi,  Eq.  (28),  and  the  standard  deviation 
of  the  phase  noise  a„  for  the  case  (i)  model.  Both  signal  and 
directional  noise  were  assumed  to  have  the  same  degree  of  random  gain 
and  phase  perturbation.  The  7 x k DOF  case  applies  to  both  hardware 
and  software  constrained  optimal  beamformers  while  the  k x k DOF 
examples  refers  only  to  the  hardware  constrained  structure. 


20 


30 


o deg 


reduction  for  case  (i). 
I (b)  a - 0.2 


B.  Interference  and  Noise  Rejection:  The  degree  to  which  noise  and 
interference  are  rejected  was  measured  by  computing  the  signal'to- 
interference  ratio  SIR  and  signal-to-total-noise  ratio  STR  observed 
at  the  beamformed  output  of  a particular  structure.  These  quantities 
were  defined  as  follows: 

A total  target  signal  output  power 

total  directional  interference  output  power 

A total  target  signal  output  power 

total  interference  and  noise  output  power 

Results  obtained  for  the  cases  studied  are  summarized  in  Fig.  8 
through  11. 

C.  Rejection  of  Target  Signal:  Figure  12  shows  the  effect  observed 

in  SIR  and  STR  when  the  mean  angle  of  arrival  of  the  target  signal 
was  other  than  bores ight.  These  curves  measure  the  penalty  incurred 
when  the  array  is  miss-pointed.  Angular  tilts  of  0 to  0.3  ir  radians 
were  investigated.  In  the  three-component  example  — i.e.  case 

(ii) — all  three  target  signals  were  assumed  to  have  the  same  offset 
angle.  It  should  be  noted  that  the  shape  of  the  conventional  beam- 
former  response  in  these  curves  is  identical  with  the  beam  pattern 
of  the  conventionally-weighted  array. 

Additional  responses  of  this  type  are  presented  in  Fig.  13  for 
case  (ii)  as  a function  of  the  amount  of  gain  noise,  in  these  plots 
the  STR  values  have  been  normalized  to  their  boresight  values  to 
better  illustrate  the  degree  of  rejection  achieved.  Similar  behavior 
was  observed  for  the  case  (i)  example. 
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Figure  8.  SIR  vs.  gain  reduction  (a)  and  phase  noise  variance  for  case  (i) 


/ 

D.  Recei ve-array  Beampatterns: 

Beam-response  patterns  were  computed  for  the  case  (ii)  example 
under  three  gain  and  phase  noise  conditions.  The  results  are  shown 
in  Fig.  14.  All  patterns  were  evaluated  at  frequency  f = 0.25  and 
the  plots  have  been  normalized  at  8 = 0 tt  radians.  As  a result, 
the  actual  mainlobe  gain  reduction  described  above  is  not  evident 
in  these  plots. 
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V DISCUSSION  AND  CONCLUSIONS 

A random  channel  model  has  been  formulated  and  presented  which 
allows  investigation  of  the  sensitivity  of  beamforming  structures  to 
random  gain  and  phase  perturbations  in  the  input  waveforms.  Procedure 
for  computing  the  appropriate  space/time  correlation  values  of  the 
received  signals  have  been  given.  These  allow  direct  computation  of 
the  sensitivity  of  any  least-mean-square  criterion  optinnal  beamformer. 
Two  specific  optimal  structures  were  studied,  each  belonging  to  the 
sidelobe  cancelling  array  class.  The  conclusions  drawn  from  this 
study  are  as  follows: 


1.  Optimal  array  processing  offers  significant  interference 
rejection  improvement  over  conventional  beamforming  even  at 
large  gain  and  phase  error  levels. 

2.  Performance  degradation  of  optimal  processors  is  largely 
dominated  by  random  phase  errors. 

3.  The  performance  of  the  hardware  constrained  optimal  beamformer 
in  the  presence  of  random  fluctuations  is  identical  to  that 

of  the  software  constrained  structure  provided  that  the  maximum 
possible  spatial  degrees  of  freedom  are  employed.  All  attempts 
to  reduce  the  number  of  spatial  degrees  of  freedom  below  this 
level  in  the  hardware  constrained  processor  led  to  poorer 
performance. 

i*.  Increasing  the  temporal  degrees  of  freedom  in  the  optimal 
processor  (i.e.,  the  number  of  taps  on  each  delay  line)  had 
no  obvious  effect  on  the  performance  of  the  processor  in 
the  presence  of  random  fluctuation. 

5.  Beam  patterns  computed  for  optimal  beamformers  show  narrower 
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widths  of  the  mainlobe  when  greater  numbers  of  spatial 
degrees  of  freedom  are  used  in  the  auxiliary  channel.  In 
addition,  these  patterns  exhibit  wider  null  depths  in  the 
presence  of  phase  noise. 
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APPENDIX  A 


To  derive  Equation  (19),  a cost  function  is  defined  as  follows: 


H(WJ 

a 


1 (w_  - w,V)'^(>  w - w \ ) , 

2 -c  -s  -a  -X-C  -s  -e'  ’ 


(A.l) 


a factor  of  1/2  is  included  to  simplify  later  arithmetic. 
Taking  the  gradient  of  (A.l), 


VH(W  ) = -W  (W  - W ) 
-a  =s  -x  c “S  -a 


(A.2) 


Letting 


VH(W  ) = 0 


W^R^W  » W ^R  W 
*s«*s«s  -a  »s  “x-c 


(A.3) 


Assuming  R^  is  positive  definite,  it  can  be  shown  that  the 
inverse  of  the  transformed  matrix  (W^R^W^^)  exists.  Then  the  optimal 


solution  of  W is. 

-a  ’ 


W opt  = (W  R W ^)‘’w  ^R  W 
-a  'as»x*s  -s  »x-c 


(A.4) 


which  is  the  resulting  optimal  formula  (19). 
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